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ABSTRACT 



Context. The detection of extrasolar planets' atmospheres requires very demanding observations. For planets that can not be spatially 
separated from their host stars, i.e. the vast majority of planets, the transiting planets are the only ones allowing to probe their 
atmospheres. This is possible from transmission spectroscopy or from measurements taken during secondary eclipse. An alternative 
is the measurement of the Rossiter-McLaughlin-effect, which is sensitive to the size of the planetary radius. Since the radius is 
wavelength-dependent due to contributions of strong planetary absorption lines, this opens a path to probe planetary atmospheres also 
with ground-based high-resolution spectroscopy. 

Aims. The major goal of our numerical simulations is to provide a reliable estimate of the amplitude of the wavelength-dependent 
Rossiter-McLaughlin-effect. 

Methods. Our numerical simulations provide phase resolved synthetic spectra modeling the partly eclipsed stellar surface during the 
transit in detail. Using these spectra we can obtain Rossiter-McLaughlin-curves for different wavelength regions and for a wavelength- 
dependent planetary radius. Curves from regions with high and low contributions of absorption lines within the planetary atmosphere 
can be compared. From these differential effects observable quantities are derived. 

Results. We applied our simulations to HD 209458. Our numerical simulations show that a detailed treatment of the limb darkening 
for the synthetic spectra is important for a precise analysis. Compared to a parameterized limb darkening law, systematic errors of 
6 msec -1 occur. The wavelength dependency of the planetary atmospheres over the NaD-doublet produce a differential effect in the 
Rossiter-McLaughlin-curve of 1.5 msec -1 for a star with a rotation velocity of 4.5 km sec -1 which increases to 4msec -1 for twice the 
rotation velocity. 

Conclusions. The Rossiter-McLaughlin-effect as tool to probe planetary atmospheres requires phase resolved, high signal-to-nose, 
high-resolution spectra taken with a stabilized spectrograph in order to obtain reliable results for slowly rotating (< 10 m sec -1 ) planet 
host stars. Stars with spectral type earlier than about F5 are a bit less demanding since the typically higher rotation velocity increase 
the amplitude of the effect to about 15 m sec -1 for a star with v sin i = 25 km sec -1 . 

Key words, planetary systems, techniques: radial velocities, Line: profiles, Stars: rotation Stars: individual: HD 209458 



1. Introduction 

Within the recent two decades, the search for extrasolar plan- 
ets has been mainly driven by the success of the radial ve- 
locity (RV) method. A large fraction of our knowledge about 
the physical properties relies, however, on the analysis of 
transiting planets. This is due to the fact that for transit- 
ing planets more information can be obtained from various 
techniques: The radius and the orbital inclination from the 
light curve analysis together with the mass function from RV- 
measurements provide the mean planetary density, an important 
constraint for the structure and evolution of the planet. The at- 
mosphere can b e probed by means of transmission spectroscopy 
during transits (|Cha rbonnea u et al.L l2002t IVi dal-Madi ar et all 
2003[l2004lBarmanll2007tlTinetti et alll2007l;lEhrenreich et all 
2008: jLecaveHer Pes Etangs etall l2008albt iPont etall 12008; 
Redfield et all |2008t Ising et all I2008a1b1: ISnellen etafl l2008t 
Swain et al.L 2008b), the albedo and thermal emission can be 
deter mined from the secondary eclipse with infrared photome- 
try (ICharbonneau et all 120051: iDeming et all 12001 120061 l2007t 
Dem orv etall l2007t iGr illmair et al.. 2007; Harrin gton et all 
2007L iRichardson et all l2007t ISwain et all l2008al) . 



From eclipsing binary systems it is known since decades, 
that the radial velocity curve shows a char a cterist i c feature dur- 
ing e clipses dSchlesingerl Il910t iRossiterl 1 19241 [McLaughlin, 
1924 . This Rossiter-McLaughlin-effect (RME) can be ex- 
plained by the partial coverage of the eclipsed star. Out of 
eclipse, the specttal lines will be symmetric since all surface el- 
ements of the (rotating) star contribute equally, thus blue and 
red shifted regions of the surface are balanced. This is changed 
during the eclipse depending on the path of the occulting body 
with respect to the eclipsed star, resulting in asymmetric lines 
due to unbalanced Doppler-shifted contributions. Depending on 
the spectral resolution of the spectrograph and the rotational ve- 
locity of the star this can be either detected as pseudo RV vari- 
ation (slow rotation or low resolution) or directly as line dis- 
tortion moving over the spectral lines during the eclipse. The 
RME has been shown to be a powerful tool for transiting plan- 
ets as well. The method has been applied for the measurement 
of the spin-orbit-alignment for eleven of the transiting planets 
dBundv & Marcvl l2000t lOueloz et all l2000t IWinn et aUj2005L 
2TM 2007; N arita etall 2007; W olf et all l2007t iBouchv et all 



2008j: ICochran et all 120081: iHebrard et all l2008t Uohnson et all 
2008tlJoshi etalll2008; Loeillet et al., 2008; Narita et al., 2008 
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With the fast growing number of bright host stars of transiting 
planets, either from ground-based wide field surveys or from 
space missions, this measurement turns into a standard tool for 
investigating properties of extrasolar planets. 

The RME was also proposed as tool for probing plan- 
etary atmospheres with ground-based observations dSnellenl 
2004). Since transmission spectroscopy as well as albedo mea- 
surements mainly rely on the high pre cision obtaina b le onl y 
from space-based measurements (see iRedfield et al.l J2008); 
ISnellen et al.l (120081) as exception), ground-based access is an 
attractive alternative. The RME therefore is a complementing 
rather than a competing method. The method detecting plane- 
tary atmospheres with the RME is based on the fact that the 
radius of an extrasolar planet depends on wavelength. This is 
caused, on the one hand, by the fact that the atmosphere of the 
planet shows variations in optical depth with respect to the wave- 
length (iBarmanl 120071: iFortnev et all 120081) . The planet will ap- 
pear slightly larger in optically thick spectral lines compared to 
wavelength regions where the atmosphere is transparent. On the 
other hand, the amplitude of the RME depends on the radius of 
the occulting body. A wavelength-resolved measurement of that 
effect provides therefore an access to the atmospheric properties 
of the transiting planet. 

While the measurements bv lSnellenl (12004 were at the limit 
of the available data, it is the scope of this paper to provide de- 
tailed modeling of the RME in order to provide a reliable anal- 
ysis tool for such demanding measurements. We therefore de- 
scribe the underlying model in Sect. 2, and present test calcu- 
lation for HD 209458 in Sect. 3, the results are summarized in 
Sect. 4. 



2. Modeling the Rossiter-McLaughlin-Effect 

Recently, there have been two papers on analytical expressions 
for the RME, expressing the RV-shift in terms of orbital parame- 
ters and radius of the planet, orbital orientation relative to the ro- 
tationa l axis, stellar rotat i on velocity, rad ius as well as limb dark- 
ening dOhta et al.. 2005; Gimenez, 2006) m aking use of earlie r 
analytical work dPetriel Il938t iHosokawal 119531: iKopall 119591) . 
While t he analytica l mode ls are well suited as fast and direct 
method, IWinn et ail (120051) argue that the the first moment of 
the spectral lines, i.e. the shift in center of gravity in wavelength 
space, as result of the analytical formula might not be identical 
to the value measured from observations. In that case, the op- 
timal match between a template and the observed spectrum is 
searched for. 

The alternative approach therefore is a finite element model 
of the star adding contributions from all surface elements to 
a synthetic spectrum which can then be treated with the same 
analys is to ols as the observat ions, e.g. as used by lOueloz et all 
(120001) and I Winn etaTI d2005l) . While computationally more de- 
manding, this approach allows to deviate from simplifications 
necessary in the analytical models allowing for a more precise 
modeling of the RME. In the following paragraphs the modeling 
of the wavelength-dependent RME will be described. 

2.1. The star 

The assumption of a spherical primary star is well suited in 
the case of a planetary secondary since deviations from spher- 
ical symmetry of the primary due to tidal interaction is negli- 
gible. Spherical symmetry also requires a slow rotation of the 
primary. For late type stars this assumption is well justified. It 



is, however, challenged when the rotation velocity becomes sig- 
nificantly larger than in the Sun. In order to prepare for such 
ca ses, we make use of the program BRUCE, originally written 
bv lTownsendl(ll997l) to model stellar pulsations in rotating stars. 
A deformation of the stellar surface is calculated as perturbation 
providing a variation of the effective temperature and surface 
gravity over the surface. As input, the stellar mass, the polar ra- 
dius and effective temperature, the equatorial velocity as well as 
a velocity amplitude due to pulsation, and the inclination i of 
the rotation axis (assumed to be constant) have to be specified. 
In our application, we set the pulsational velocity to zero (this 
could, however, be included in order to investigate a contribution 
of stellar pulsation to RV noise). In an output file, BRUCE pro- 
vides effective temperature, surface gravity and projected radial 
velocity as well as the cosine of the angle between the normal 
direction of the surface element and the line of sight of the ob- 
server (fj.) for the required number of surface elements (typically 
50000 in our test calculations). 

2.2. Synthetic spectrum of the host star 

In previous modeling authors used template spectra, e.g. the so- 
lar spectrum, or mean intensities from model atmosphe res as 
input spectra. In b oth cases, a limb darkening law (e.g. IClaretl 
d2000ll2003Ll2004l) ) has to be applied. For a detailed wavelength- 
dependent modeling of the RME the available broad-band limb 
darkening coefficients are not sufficient, instead this would re- 
quire determination of limb darkening coefficients from model 
atmospheres adapted to the spectral resolution of the observa- 
tions. In order to avoid errors due to the parameterization of the 
limb darkening law, specific intensities for all surface elements 
of the star can be used directly for a synthetic spectrum. 

The input spectra were synthesised using specific intensities 
(energy emitted in a given direction n per time, unit area, fre- 
quency and solid angle interval) from a grid of PHOENIX mod- 
els. The grid covers the main sequence in steps of 200 K for 
the effective temperature and with two surface gravities (log g=4 
and log g=4.5). The models have been calculated for spherically 
symmetric atmosph eres under the assumption of local thermody- 
namic equilibrium (Hauschildt & Baron, 1999; Hausc hildt et all 
1999). The general setup is described in iBrott & Hauschildt! 
(2005 1), however we are adopting the revised solar abundances 
of lAsplund et a l. (2005) to calculate the chemical composition 
(for solar metallicity in the present case, but being available for 
other metallicities as well). Specific intensity spectra for each 
model computed at a resolution of 0.01 A for different values of 
the angle 6 between the normal direction of the atmosphere and 
the direction of the line of sight (p. = cos 9). 19 intensities per 
wavelength point are extracted from all core rays and a subset 
of the tangent rays to sample the range < fi < 1 evenly. Thus 
while our model treats the global geometry of the stellar surface 
exactly including rotational flattening, it is locally approximated 
in spherical geometry, while interpolating to the exact effective 
temperature and gravitational acceleration for each atmospheric 
segment. Since the deformation and the total curvature are very 
small for dwarf stars, these approximations lead to negligible 
errors. The spectrum is then synthesized by linear interpolation 
within the available grid of intensities on a logarithmical equidis- 
tant wavelength grid of the required spectral resolution and shift- 
ing to the required Doppler velocity. For comparison, we can 
also calculate intensities from the disk-averaged spectrum by ap- 
plying a linear limb darkening law. The undisturbed spectrum of 
the planet host star is now obtained and can be used as template 
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for cross-correlations to determine the RME and it serves as ref- 
erence spectrum of the next step. 



2.3. The Planet 



The RME is dependent on the size of the planet as well 
as on its orbital parameters. Since transiting planets are 
typically close in, the assumption of a circular orbit is 
well justified, however, eccentric orbits with eccentricities 
above 0.1 have b een detected i n four transiting planets: 
GJ436b, ecc=0 1 5 iDeming et al.1 (120071) . XO-3 b, ecc=0.26 
Uohns-Krull etaD d2008l), HP 14 7506b, ecc=05163lBakos et alJ 
d2007|), HD 17156b, ecc=0.6753 lFischer et alJ d2007|). We there- 
fore specify the orbital separation and period. An extension to- 
wards eccentric orbits would be possible in case this is needed. 
The orientation of the planetary orbit with respect to the line 
of sight towards the observer also has to be specified since the 
Doppler shifts due to stellar rotation breaks the symmetry in this 
case. The orientation of the orbital rotation axis relative to the 
line of sight needs additional two parameters, the inclination i 
and the tilt of the orbit in the plane perpendicular to the line of 
sight A. 

Synthetic spectra can now be calculated. We specify an 
equidistant phase sampling (starting shortly before the first con- 
tact and ending shortly after the fourth contact), the radius of the 
planet as well as the angles i and A. For each orbital phase, the 
elements on the stellar surface which are eclipsed are identified 
and their spectra are interpolated as in the previous step. Their 
contribution is then subtracted from the reference spectrum re- 
sulting in the synthetic spectrum for the given orbital phase. This 
procedure can be repeated for different planetary parameters. 




Fig. 1. Wavelength-dependent planetary radius of HD 209458b 
dBarmanl |2007|) (dashed line) and smoothed over 100 A bins 
(full line). 




2.4. Wavelength-dependent planetary radius 



Fig. 2. Rossiter-McLaughlin-Effect for HD 209458 (orbital mo- 
tion subtracted): Top panel: Planetary radius Rf>=1.37Rj up and 
Rf=1.31Rj up (thin lines, full simulation), Rp=1.34Rj up (thick 
line, full simulation) and Rp=1.34Rj up (++++, interpolation 
procedure using Rp=1.37 Rj up and Rp=1.31 Rj up ). Lower panel: 
Difference between the full simulation and interpolation proce- 
dure for Rf=1.34Rj up . See text for details. 



From radiative transport calculations, Barman (2007) deduce a 
wavelength dependency of the planetary radius as the result of 
the varying optical depth of the planetary atmosphere (see FigQ] 
for an adopted version of his Fig.l). From his results a variation 
in the order of 5% can be expected, with a promi nent peak e.g. 
at the resonance lines of sodium. Predictions of iFortnev et al.l 
(2008) indicate a variation between 2% and 5% depending on the 
class of irradiated atmosphere. In order to take the wavelength- 
dependent radius into account, we calculate two sets of synthetic 
spectra (F up , F\ ow ), one with the upper limit and one with the 
lower limit of the radius (r up , ri ow ). For each wavelength point, 
we then interpolate the two synthetic spectra for the giv en wave- 
length point to the value according to the prediction of iBarmanl 
d2007l) . i.e. F inU = F iovl + (r inU - r iovj ) ■ (F up - F iovj )/(r uv - r iov/ ). 
It should be noted, that we are not restricted to these prediction 
but can accommodate any given wavelength-depend ent radius, 
e.g. b ased on observations of transmission spectra Sing et al. 
(2008a). The application of a linear interpolation is well justi- 
fied as can be seen in Fig. [2] In the upper panel we display the 
RME for the lower and upper limit of the predicted wavelength- 
dependent radius as well as the RME for an intermediate value, 
one calculated directly as described in the previous step the other 
one interpolated between the two limiting cases. The lower panel 
shows the difference between the two, which are of the order of 
cm sec -1 and therefore certainly below current detection limits. 



3. Results: Test case HD 209458 

In the following, we present test calculations of the RME 
and transit light curves for HD 20 9458. The stellar and p lane- 
tary parameters are adopted from iNordstrom et al.l d2004l) and 
iBrown et al.l d200ll) . namely we use R» iPO iar = 1.15 R , M* = 
1.05 M , vsin/ = 4.5kmsec~ I , T e ff, po iar=6000K, and an incli- 
nation of the rotational axis of 90° as well as an orbital separa- 
tion of 0.047 AU, an orbital inclination of 86.6° and no tilt of the 
orbit in the plane perpendicular to the line of sight. The plane- 
tary radius is treated wavelength-dependent (FigJTJ as d escribe d 
in the previous section using the predictions of Barman (2007). 

3.1. The effect of limb-darkening 

The form of the RME is not only dependent on the stellar and 
planetary parameters, it is also sensitive to the treatment of limb 
darkening. As described in Sect. 2, our spectral synthesis re- 
lies on the co-addition of specific intensities for each surface 
element, interpolated to the actual value of the angle between 
the normal of the element and the line of sight. In Fig. [3] and [4] 
we compare the averaged RME from a reference region (5000- 
5200 A) from our full simulation including specific intensities 
to two versions of simulations using the flux and a linear limb 
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Fig. 3. The Effect of limb darkening: Top panel: Rossiter- 
McLaughlin-Effect (orbital motion subtracted) for the wave- 
length region 5000-5200 A from synthetic spectra using specific 
intensities (full line) compared to the application of a linear limb 
darkening law (dashed line). Lower panel: Difference between 
the full simulation and the limb darkening approximation. See 
text for details. 
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Fig. 4. As Fig. [3] but for the application of a wavelength depen- 
dent linear limb darkening law shown in Fig. [5] 



darkening law I^j = 1 + (// - 1)m,i for the spectrum synthesis. In 
th e first case we u se a constant coefficient of u A — u — 0.58 used 
bv lSnellenld2004l) . in the second case we use a polynomial fit of 
second order to the wavelength-dependent linear limb darkening 
coefficient (Pig.[5j>. 

The deviation between the simplified limb darkening ap- 
proach and the detailed simulation using specific intensities is 
most evident when the planet covers the stellar limb, i.e. before 
the second and after the third contact. Here, the differences be- 
tween the limb darkening law and the actual /i-dependence of the 
specific intensities are largest. A limb darkening law therefore 
mocks a star with a different effective radius. The differences can 
be reduced with a wavelength-dependent limb darkening (Figj4j 
law but not completely suppressed. 

Systematic differences between the detailed treatment and 
a limb darkening approximation are up to 6 msec 1 in the first 
case and 2 m sec -1 in the second case. A detailed treatment of the 
RME might therefore help to reduce the error bars in determina- 
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Fig. 5. Wavelength-dependent linear limb darkening coefficient 
determined from the applied PHOENIX models smoothed over 
100 A bins (full line) and polynomial fit of second order (dashed 
line). 



tions of the spin-orb it alignments or misalignements of transiting 
extrasolar planets dBundv & Marcvi l2000t lOueloz et all [2 000: 
I Winn et all I2005L 120061 120071 iNarita et all 120071: IWolf et al. 
2007 1 : iBouchv et all 120081: ICochran et all 120081: iHebrard et al. 
2008; iJohnson et all 120081: iJoshietall 120081: lLoeillet et al' 
2008; Na ritaet alll2008l) . Due to the wavelength dependency, a 
detailed treatment seems also necessary in order to avoid confu- 
sion between effects of a wavelength-dependent planetary radius 
and a wavelength-dependent limb darkening. 



3.2. The Wavelength-dependent Rossiter-McLaughlin-Effect 

With our simulations and predictio n s for the wavelength depen- 
dency of planetary radii (Barmarjl 120071: iFortnev et all [2008) 
it is now possible to simulate the wavelength dependency of 
the RME and to estimate the observability. All test calcula- 
tions were applied to HD 209458. It should, however, be noted 
that our simulations are not restricted to that object. We as- 
sume a wavelength range from 5000-6500 A, roughly covered 
by a spectrograph stabilized by an iodine-cell and covering the 
NaD-doublet where the planetary radius is increased by about 
5% due to absorption within the planetary atmosphere. We used 
the phase-dependent synthetic spectra as described in Sect. 2, 
performed cross-correlations stepping through the spectral re- 
gion with a step size of 25 A, and a width of 100 A, i.e. averag- 
ing the RME over each 100 A-bins. The correspondingly aver- 
aged wavelength-dependent radius of HD 209458b can be seen 
in Fig. 03 As a reference, we use the RME averaged over the 
region 5000-5200 A, because the predicted radius is least influ- 
enced by absorption in the planetary atmosphere in our spectral 
region of interest. The differential effect is illustrated in Fig. [6] 
where we show the RME of the reference region, the 100 A-bin 
centered around NaD, as well as the difference between the two. 
The amplitude of the effect is about 1 .5 m sec 1 , about the quan- 
tity derived by Snellen (2004). It should be noted that the differ- 
ential effect is qualitatively different from that due to simplifica- 
tions using parameterized limb darkening laws (Figs.[3]and|4]i. 
The maximum effect occurs when the planet covers most of the 
blue- or red-shifted hemisphere, i.e. after the second and before 
the third contact. The wavelength dependency of the planetary 
radius therefore has a different signature in the RME than limb 
darkening. Using a simplified limb darkening model as compari- 
son for the observations would in principle allow to discriminate 
between effects of the limb darkening and the wavelength de- 
pendency, in practice the necessary accuracy would be difficult 
to achieve. A detailed model is therefore necessary to obtain re- 
liable results. 
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Fig. 6. Rossiter-McLaughlin-Effect (orbital motion subtracted) 
for the wavelength region 5840-5940A (100A bin centered 
around the NaD-doublet): Simulation with the wavelength- 
dependent radius for that wavelength interval (top panel, full 
line) versus a fixed radius Rp = 1.31 Rj up (top panel, dashed 
line) and difference between the two (lower panel). 



In Fig. [7] we show the wavelength dependency of these dif- 
ferences: The amplitude of the difference between the RME of 
each 100 A-bin and the RME of the reference region is plotted 
against the central wavelength of each bin as full line. While this 
indicates an signal of about 3msec~' with a peak at the NaD 
doublet, it should be noted that about half of the effect is due to 
the wavelength dependency of the limb darkening. This can be 
seen from the dotted line, where we keep the planetary radius 
constant at Rp = 1.31 Rj up , the wavelength-dependent size cor- 
responding to the reference region. Here we can also see a peak 
at NaD, howev er, with a low er amplitude (1.5 m sec" 1 ). The ob- 
servations of ISnellerj d2004l) indicate a measured amplitude of 
1.7 msec -1 for HD 209458 from UVES spectra, in accordance 
to our simulations. 

As a possible strategy for an analysis we therefore suggest 
to use the observed RME of the reference region to obtain the 
necessary stellar and planetary parameters from a fit of the sim- 
ulated RME (stellar radius, orientation of the orbit with respect 
to the rotation axis of the star as well as with respect to the line 
of sight and especially the planetary radius at this wavelengths) 
and use the corresponding synthetic RME as reference. The dif- 
ferential effect as displayed in the lower panels of Fig.|6]and|7] 
would then be the result of the analysis. 

3.3. The influence of the stellar rotation velocity 

For a broader application of these technique, the influence of the 
rotation velocity of the host star is the most critical parameter. 
While HD 209458 is a slow rotator, several of the known host 
stars of transiting planets rotate faster. In order to avoid confu- 
sion between effects of changing stellar parameters for different 
host stars, we compared the results from the previous paragraph 
now to simulation for models where all parameters are kept fixed 
except the rotation velocity. The results are presented in Fig. [8] 
As expected, the amplitude of the observable effect increases 
with the rotation velocity of the host star. While the amplitude 
is about 1.5msec _1 for a rotation velocity of 4.5kmsec~ 1 it in- 
creases to about 4msec~ 1 for a rotation velocity of 9 km sec -1 . 
For planet host stars with spectral types earlier than about F5 ro- 



Fig. 7. Wavelength-dependent Rossiter-McLaughlin-Effect at 
the phase of maximum amplitude (-0.012): Top panel: 
Difference between the wavelength-dependent RME and 
the RME from the reference region 5000-5200 AFull line: 
Simulation with the wavelength-dependent radius, dotted line: 
Simulations with a a fixed radius Rp = 1.31 Rj llp . Lowe panel: 
The difference between the two simulations from the top panel. 
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Fig. 8. The effect of stellar rotation: Same as lower panel 
of Fig. [7] but for stellar rotation velocities of v eq = 
4.5, 6.0, 9.0, 25.0km sec -1 (thin, dotted, dashed, and thick line). 



tation is significantly larger compared to those of later spectral 
type. Bright planet host stars of earlier spectral types are there- 
fore the primary targets for this method. 

3.4. Wavelength-dependent transit light curves 

As a by-product, we can also simulate transit light curves by in- 
tegrating the flux for each phase point. We can therefore also 
provide detailed simulations for transmission photometry. As an 
example, we show in Fig.[9]the transit light curve for HD 209458 
for the reference region 5000-5200 A (assuming a box-like fil- 
ter) and for the average of two 3 A-bins each centered around 
one of the NaD components. The difference between the two is 
provided in the lower panel as possible observational result, sim - 
ilar to those o fjCharbon neau et alJd2002l) : lRedfield et alj ([2008); 
ISnellen et al] ((20081) : ISing et all d2008ai) . The differential effect 
is largest during ingress and egress because the slightly larger 
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Fig. 9. Wavelength-dependent transit light curves: average from 
two wavelength regions centered around the two NaD doublet 
components of 3 A width each (top panel, full line), from a refer- 
ence region 5000-5200 A (dashed line), and difference between 
the two (lower panel). 

planetary radius at the NaD-lines causes a slightly early first 
and slightly later fourth contact. The differential effect during 
phases between the second and third contact is a combination of 
the slightly deeper transit due to a larger radius at NaD and the 
differential effect of the limb darkening between the two wave- 
length regions. This part of the curve therefore depends critically 
on the wavelengths of the reference region. 

Since (ground-based) measurements for the differential ef- 
fect of the transit light curves with wavelength regions of 
strong absorption lines of the planetary atmosphere rely on 
high resolution and high-quality spectra, both, the RME and the 
wavelength-dependent transit light curve, can be evaluated si- 
multaneously. With our simulations as comparison, the detection 
efficiency can therefore be increased significantly. 

4. Conclusion 

We performed simulations of the Rossiter-McLaughlin-effect 
using a predicted wavelength dependency for the transiting 
planet in order to evaluate the possibility to detect atmospheres 
of extrasolar planets with ground-based high-resolution spec- 
troscopy. We showed that a detailed treatment of the spectrum 
synthesis of the partly eclipsed stellar surface is important in or- 
der to detect these subtle effects. In typical applications for the 
measurement of spin-orbit alignment, we obtain systematic dif- 
ferences up to 6 ms 1 in simulations tuned for HD 209458. A 
faster rotation of the primary star would increase these system- 
atic effects. 

The slightly larger planetary radius at the NaD-resonance 
line, caused by additional absorption within the planetary at- 
mosphere, produces a slightly larger RME compared to a wave- 
length region where the planetary atmosphere hast a very small 
contribution (e.g. the region 5000-5200 A). Our simulations for 
HD 209458, predict an amplitude of the effect of 1 .5 m s _1 . This 
is a very small effect which requires several phase-resolved high 
resolution spectra from several transits taken with a stabilized 
spectrograph. While the rotation velocity of HD 209458 is low 
(4.5 km s" 1 ), systems with a primary of spectral type F5 or earlier 
typically show higher rotation velocities. For rotation velocities 
of 25 km s _1 the amplitude of the effect is already 15 m s -1 . With 



increasing rotational velocity, the accuracy of the radial velocity 
measurements will, however, be reduced due to the broadening 
of the stellar lines. 

From these simulation, we can also derive wavelength- 
dependent transit light curves. A larger planetary radius at the 
NaD-resonance line could be detected from a comparison be- 
tween the light curve at this wavelength to one at a reference 
region with small contributions from the planetary atmosphere. 

During the recent years, atmospheres of extrasolar planets 
have been probed with different techniques. Photometric ob- 
servations of the planetary albedo and secondary eclipses has 
to rely on space-based IR instruments, while the transit spec- 
troscopy has recently been successful also from the ground (see 
Sect. 1). In this paper it has been been shown that the analy- 
sis of the wavelength-dependent Rossiter-McLaughlin-effect is 
an alternative approach for a ground-based detection. While ac- 
cess from the ground will potentially enable a broader applica- 
tion, the information of the ground based transit spectroscopy 
as well as our proposed application of the RME is restricted to 
strong absorption lines in the planetary atmosphere. An advan- 
tage of the wavelength-dependent RME is that it can be obtained 
as a by-product of either ground-based transit spectroscopy, pro- 
viding an independent detection of absorption lines in plane- 
tary atmospheres, or as by-product to RME measurements for 
spin-orbit alignment provided that the requirement for the wave- 
length stability of the spectrograph is high enough to measure 
the wavelength-dependent effect. 
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